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y  We  show  that  under  realistic  ionospheric  conditions,  barium  ion  clouds  hundreds  of  meters  in  diameter 
can  be  long-lived,  quasi-stable,  nonbifurcating  structures.  These  structures  may  resemble  “tadpoles,”  with 
a  dense  head,  steep  density  gradients  at  the  front,  and  a  long,  less  dense  tail.  We  assume  that  these  struc¬ 
tures  are  the  final  products  of  the  recursive  bifurcation  of  a  considerably  larger  barium  ion  cloud,  i.e., 
striations.  The  realistic  ionospheric  conditions  to  which  we  refer  consist  of  a  barium  ion  cloud  with  ion 
temperatures  Tj  of  approximately  1000°  K,  coupled  electrically  to  a  background  ionosphere  of  lower 
compressibility  than  itself,  i.e.,  and  F  region.  We  show  analytically  that  this  combination  of  finite  Tj  and 
relatively  incompressible  background  results  in  an  effective  diffusion  of  barium  plasma,  but  more 
importantly,  of  total  magnetic-field-line-integrated  Pedersen  conductivity  ,'S^C.  The  diffusion  coefficient  has  a 
special  form  which  allows  the  inner  portions  of  the  cloud  to  diffuse  slowly,  giving  the  cloud  a  long  lifetime, 


;q  D.S’RiBu' 

n  a  w  a  .  i 

JNCI  ASS'  F  •*  <“> 

N  L  Ml  ' 

12*  NAVE  }F  AeSA  U'vlS'.A. 


(Continues) 


.'I  AdS'-AA.ZT  "  ZEASS.F  CATION 

UNCLASSIFIED 


;i:  'F  JF-'  NF  Ni.VflfcFi  22c  ZFL'Cl  S.VIIGl 

In,  mu.-  »  ii-  .ij.- 

(202)  767-3215 


S.  T.  Zaleask 


Code  4780 


11.  TITLE  (Include  Security  Classification) 


Finite  Temperature  Effects  on  the  Evolution  of  Ionospheric  Barium  Clouds  in  the  Presence  of  a 
Conducting  Background  Ionosphere 

I.  A  High  Altitude  Incompressible  Background  Ionosphere 
16.  SUPPLEMENTARY  NOTATION  (Continued) 

This  research  was  sponsored  by  the  Defense  Nuclear  Agency  under  Subtask  S99QMXBC, 
work  unit  00102  and  work  unit  title  “Plasma  Structure  Evolution.” 


19.  ABSTRACT  (Continued) 

and  allows  the  outer,  less  dense  portions  of  the  cloud  to  diffuse  rapidly,  preventing  cloud  bifurcation. 
Numerical  simulations  of  the  full  nonlinear  dynamics  are  then  used  to  show  that  this  diffusion  does  in 
fact  give  rise  to  quasi-stable  barium  striations  hundreds  of  meters  in  diameter.  These  findings  are 
consistent  with  the  linear  analysis  of  Francis  and  Perkins  (1975). 
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FINITE  TEMPERATURE  EFFECTS  ON  THE  EVOLUTION  OF  IONOSPHERIC  BARIUM 
CLOUDS  IN  THE  PRESENCE  OF  A  CONDUCTING  BACKGROUND  IONOSPHERE 

I.  A  High  Altitude  Incompressible  Background  Ionosphere 
1.  Introduction 

The  nonlinear  evolution  of  an  artificial  plasma  cloud  (e.g.,  a  barium 
ion  cloud)  released  in  the  earth's  ionosphere  is  characterized  by  an 
overall  bulk  JS  x  motion,  a  one-sided  steepening  of  plasma  gradients,  and 
by  bifurcation  [see  Ossakow,  1979;  Ossakow  et  al.,  1982  and  references 
therein]  a  process  which  we  define  here  to  mean  the  splitting  of  the  cloud 

into  two  or^  more  pieces.  These  pieces  are  sometimes  referred  to  as 

"striations".  All  of  the  above  motion  and  the  resultant  plasma  structures 

are  observed  to  be  magnetic-field-aligned.  That  is,  the  plasma  motion 
parallel  to  magnetic  field  lines  consists  primarily  of  diffusion  and 
falling,  while  the  perpendicular  plasma  motion  is  virtually  identical  for 
all  plasma  along  a  given  field  line.  Hence  one's  best  view  of  the 
structuring  of  the  cloud  is  from  a  vantage  point  looking  parallel  to  the 
magnetic  field  jJ.  This  simplified  description  of  barium  cloud  dynamics, 
seen  along  a  line  of  sight  parallel  to  B_,  is  depicted  in  Fig.  1.  The 

sequence  of  drifting,  steepening  and  bifurcating  i3,  for  sufficiently  large 
striations,  recursive:  each  striation  becomes  a  new  cloud  which  in  turn 
forms  its  own  striations.  In  the  absence  of  some  dissipative  mechanism,  it 
is  believed  that  this  process  would  continue  indefinitely.  In  reality,  the 
striation  size  becomes  small  enough  to  be  comparable  to  the  characteristic 
scale  lengths  of  physical  dissipation  mechanisms,  and  further  bifurcation 
ceases.  A  previous  study  by  McDonald  et  al.  [1981]  attempted  to  quantify 
this  concept  within  the  context  of  a  "one-level  model",  a  mathematical 
model  in  which  all  conducting  plasma  was  constrained  to  lie  in  a  single 

two-dimensional  plane  perpendicular  to  the  magnetic  field.  Implicit  in 
Manuscript  approved  March  15,  1984. 
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this  model  is  the  assignment  of  a  single  Pedersen  mobility  and  hence 
compressibility  to  both  the  barium  and  background  ionospheric  plasma.  In 
this  model*  the  leading  candidate  for  a  dissipative  mechanism  is  the 
ambipolar  diffusion  of  electrons  caused  by  electron-ion  collisions.  Using 
this  mechanism,  McDonald  et  al.  [1981)  calculated  final  unbifurcating 
strlatlon  sizes  (from  a  so-called  U  shaped  curve)  for  barium  clouds  of  ~  20 
m  (perpendicular  to  _B).  However,  there  is  experimental  evidence  for  the 
existence  of  long-lived  "blobs"  of  plasma  several  hundred  meters  in 

diameter.  The  apparent  stability  of  such  large  plasma  structures  has  yet 
to  find  a  satisfactory  explanation  within  the  one-level  model.  Ossakow  et 
al.  [1981]  proposed  including  a  zero  temperature  second  level  compressible 
background  ionosphere  below  the  cloud  (see  Scannapieco  et  al.,  1976).  This 
would  allow  image  striations  to  build  up  and  allow  the  conductivity  in  a 
strlatlon  to  be  amplified.  This  In  turn  would  allow  for  larger 

conductivity  ratios  than  if  one  had  just  a  one  level  cloud.  In  turn,  this 
could  result  in  kilometer  size  unbifurcating  striations  by  extrapolating 
the  U  shaped  curve  of  McDonald  et  al.  [1981]  to  higher  conductivity  ratios. 

Here  we  postulate  that  the  existence  of  a  conducting  background  plasma 
with  ion  Pedersen  mobility  different  from  that  of  the  barium  plasma, 
combined  with  the  finite  temperature  of  the  barium  plasma,  may  allow  the 
existence  of  stable  striations  with  scale  sizes  of  hundreds  of  meters.  The 
essence  of  the  mechanism  is  a  phenomenon  commonly  referred  to  as  "end 
shorting"  first  investigated  in  the  context  of  barium  clouds  by 
Shiau  and  Simon  [1974],  who  considered  a  completely  incompressible 
background  ionosphere,  and  showed  that  the  normal  electron  ambipolar 
diffusion  rates  would  be  replaced  by  some  fraction  of  the  much  higher  ion 

diffusion  rates  caused  by  collisions  of  Ions  with  neutrals.  Francis  and 

Perkins  [1975]  generalized  the  work  of  Shiau  and  Simon  [1974]  by 


considering  the  cases  of  both  incompressible  and  compressible  background 


ionospheres.  They  concurred  with  Shiau  and  Simon  [1974]  that 

incompressible  backgrounds  exert  a  stabilizing  influence,  but  noted  that 
compressible  backgrounds  are  destabilizing .  In  contrast,  Vickrey  and 
Kelley  [1982]  concluded  that  a  background  (compressible)  E  region  should 
exert  a  stabilizing  influence  on  high  latitude  irregularities.  We  shall 
address  the  question  of  the  correctness  of  the  above  two  analyses  in  a 
later  paper.  For  the  moment  we  note  that:  1)  Francis  and  Perkins  [1975] 
performed  a  rigorous  stability  analysis  of  the  problem;  Vickrey  and  Kelley 
[1982]  did  not;  2)  Francis  and  Perkins  [1975]  self-consistently  included 
image  formation  (compressibility)  in  both  the  cloud  and  the  background; 
Vickrey  and  Kelley  [1982]  did  not  (they  took  a  passive  load  model);  3) 
Vickrey  and  Kelley  [1982]  included  recombination  chemistry;  Francis  and 
Perkins  [1975]  did  not. 

Our  own  conclusions  are  based  on  our  view  of  barium  cloud  dynamics  as 
consisting  of  a  two-stage  feed-back  loop: 

1)  At  any  given  time,  the  distribution  of  plasma  density  will, 
through  the  effect  of  these  distributions  on  magnetlc-fleld-llne  integrated 
Pedersen  and  Hall  conductivities,  bring  about  the  creation  of  polarization 
electric  fields  whose  purpose  it  is  to  maintain  quasi-neutrality; 

2)  these  electric  fields,  through  Hall  and  Pedersen  mobilities, 
affect  the  velocity  of  the  plasma,  which  in  turn  affects  the  distribution 
of  plasma  density  at  the  next  instant  of  time. 

We  emphasize  that  it  is  only  through  its  influence  on  magnetic-field- 
line  integrated  Pedersen  and  Hall  conductivities  that  a  change  in  plasma 
distribution  [e.g.,  diffusion]  will  affect  the  polarization  fields,  which 
are  the  engine  of  plasma  structure. 


In  the  present  study  we  have  isolated  the  "enhanced”  difusive  effects 
(due  to  finite  temperature  effects  and  a  two  level  model,  where  only  ion- 
neutral  collisions  are  included)  on  barium  cloud  striation  evolution.  A 
major  conclusion  derived  from  the  nonlinear  numerical  simulation  results 
presented  here  is  that  coupling  to  a  nearly  incompressible  background 
ionosphere,  i.e.,  an  F  region,  can  result  in  a  cessation  of  striation 
bifurcation  in  F  region  ionospheric  barium  clouds. 

In  addition,  a  simplified  analysis  of  the  equation  for  the  time 
evolution  of  the  total  field  line  integrated  Pedersen  conductivity  shows 
that  coupling  to  a  compressible  background  ionosphere,  i.e.,  E  region, 
would  result  in  destabilization  of  the  striations  (see  Eq.  (4-36)  and  the 
discussion  following).  Both  of  these  results  are  consistent  with  the 
linear  analysis  of  Francis  and  Perkina  [1975]. 

In  section  2  we  describe  the  motion  of  ionospheric  plasma  and  In 
section  3  the  mathematical  model  is  presented.  Section  4  presents  a 
simplified  two  level  model.  Section  5  describes  the  diffusion 
characteristics  for  a  simplified  two  level  model  with  a  nearly 
incompressible  background.  Numerical  simulation  results  using  an 
incompressible  background  are  presented  in  section  6.  Section  7  presents 
the  conclusions  and  thoughts  on  future  work. 


2.  The  Motion  of  Ionospheric  Plasma 

We  shall  be  concerned  here  with  the  motion  of  plasma  consisting  of 
ions  and  electrons  in  the  presence  of  a  neutral  gas  and  magnetic  field  B, 
subject  to  an  external  force.  We  shall  also  be  interested  in  the  electric 
current  arising  from  the  differential  motion  of  the  various  species 
comprising  the  plasma.  In  the  course  of  deriving  the  equations  we  make  the 
following  assumptions: 

1)  We  assume  the  plasma  can  be  adequately  described  by  the  fluid 
approximation.  This  assumes  that  the  effective  collision  rate  of  each 
plasma  species  with  itself  is  sufficiently  high  to  maintain  near  Maxwellian 
distribution  functions  on  time  scales  short  compared  to  the  times  of 
interest,  and  is  well  satisfied  for  the  plasmas  we  treat  here. 

2)  We  assume  that  the  electric  fields  E_ are  electrostatic  (i.e., 

V  x  E  *  0)  and  hence  can  be  described  using  a  scalar  potential  $  such  that 
_E  ■  -  V<b .  Note  that  this  Implies  3B/3t  -  0.  The  validity  of  this 
assumption  can  be  related  to  the  fact  that  the  Alfven  velocity  is  much 
larger  than  any  other  propagation  speed  of  interest  for  the  plasmas  we 
treat  here.  The  assumption  Is  also  checked  a  posteriori  by  verifying  that 
the  calculated  currents  and  displacement  currents  produce  negligible  time 
variations  in  which  in  turn  produce  negligible  V  x  E. 

3)  We  assume  plasma  quasi-neutrality;  that  is, 

l  n1<l1  "  nee  (2.D 

where  n  is  the  number  density,  q  is  ion  species  charge,  e  is  the  electron 
charge,  the  subscripts  i  and  e  refer  to  ions  and  electrons  respectively, 
and  the  sum  13  taken  over  all  ion  species.  This  assumption  is  a  statement 
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that  the  Debye  length  is  small  compared  to  all  length  scales  of  interest, 
and  again  can  be  verified  a  posteriori  by  evaluating  V  •  E.  Note  that  this 
assumption  implies  that  7  •  J_  =  0,  where  J  is  the  electric  current. 

In  addition  to  the  above  there  are  some  other  assumptions  which,  while 
they  are  not  essential  to  the  basic  model,  are  nonetheless  valid  for  many 
of  the  physical  situations  which  we  shall  treat  and  impart  a  simplicity 
which  we  shall  find  convenient  here. 

4)  We  assume  the  electrostatic  potential  $  to  be  constant  along 
magnetic  field  lines.  As  we  shall  see  later,  the  electrical  conductivity 
along  magnetic  field  lines  is  much  greater  than  that  perpendicular  to 
magnetic  field  lines,  meaning  that  appreciable  differences  in  potential 
along  a  field  line  will  quickly  be  reduced  by  the  resultant  current.  This 
assumption  will  break  down  for  sufficiently  small  scale  lengths 
perpendicular  to  the  magnetic  field,  and  for  sufficiently  large  distances 
along  the  magnetic  field. 

5)  We  assume  that  the  inertial  terms  in  the  plasma  species  momentum 
equations,  i.e.,  the  left  hand  side  of  Equation  (2.3),  are  negligible  with 
respect  to  the  other  terms  in  the  equation.  This  assumption  is  justified 
whenever  the  time  scales  of  interest  are  longer  than  the  mean  time  between 
collisions  for  ions. 

6)  We  neglect  collisions  between  ions  and  electrons  and  between  ions 

of  different  species.  Later  we  shall  also  neglect  collisions  of  electrons 

with  neutrals.  There  are  two  reasons  for  this.  First,  the  ion-neutral 

collision  term  v  can  be  shown  to  be  the  dominant  term  in  the  diffusion 

m 

physics  we  consider  here.  (One  must  be  careful,  however,  since  these  same 
assumptions  will  yield  zero  diffusion  when  the  background  conductivity  is 


set  to  zero.  To  obtain  the  correct  electron  ambipolar  rate  in  this  limit, 


the  electron-ion  and  electron-neutral  collision  terms  must  be  retained.  We 


are  not  interested  in  this  limit  in  this  paper).  Second,  although  exact 
closed  form  expressions  for  the  ion  and  electron  velocities  in  terms  of  the 
applied  forces  is  possible  when  ion-electron  collisions  are  retained 
[ Fedder ,  1980] ,  the  expressions  are  considerably  more  complex  than  those 
given  below. 

The  continuity  and  momentum  equations  describing  a  sir  e  ion  species 
and  its  associated  electrons  are: 


3n 

~  +  V  •  (n  v^)  =  0 


(2.2) 
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where  the  subscript  a  takes  on  values  i  and  e  (denoting  ions  and  electrons, 
respectively),  n  is  the  number  density,  _v  is  the  fluid  velocity,  P  is 
pressure,  E_  is  the  electric  field,  is  the  gravitational  acceleration,  q 
is  the  charge,  \>an  is  the  collision  frequency  with  the  neutral  gas,  is 
the  neutral  wind  velocity,  c  is  the  speed  of  light,  and  m  is  the  particle 
mass.  We  can  rewrite  this  equation  as 


F  /m  + -  (v  x  B)  -  v  v  = 

—a  a  me  — a  —  an  —a 
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If  we  place  ourselves  In  a  Cartesian  coordinate  system  in  which  _B  is 

aligned  along  the  z  axis,  and  if  we  treat  as  a  given  quantity  then  a 

componentwise  evaluation  of  Equation  (2.4)  yields  a  set  of  three  equations 

in  three  unknowns,  the  three  components  of  v  .  The  formal  solution  is 
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The  vector  subscripts  1  and  ’’  refer  to  the  components  of  the  vector  which 


are  perpendicular  and  parallel  respectively  to  z .•  The  quantities  k.  , 


and  kQ  above  are  referred  to  as  the  Pedersen,  Hall,  and  direct  mobilities 

respectively.  It  should  be  pointed  out  that  Equations  (2.6)  and  (2.7)  are 

only  truly  closed  form  expressions  when  the  inertial  terms  (the  last  term 

on  the  right  hand  side  of  Equation  (2.5))  are  neglected,  an  assumption  we 

have  made  previously.  Typical  ranges  for  collision  frequencies 

are:  V£n  ~  30  sec-*-,  ven  ~  800  sec~l  at  150  km  altitude;  and  v^n  ~ 

10  ^  sec  \  v  ~  1  sec  ^  at  500  km  altitude, 
en 

As  we  will  see  later,  we  will  use  the  concept  of  "layers"  to 
distinguish  the  various  ion  species,  so  for  the  moment  we  consider  only  a 
single  ion  species,  denoted  by  subscript  i,  and  the  associated  electrons, 
denoted  by  subscript  e.  We  also  consider  only  singly  charged  ions  so  that 
q^  *  e  and  q  -  -  e.  Noting  that  ven/fie  *  0,  we  obtain 
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We  now  define  the  perpendicular  current 


-1  -  L  a 
a 

Substituting  Equations  (2.13)  through  (2.15)  and  (2.6)  into  Equation 
(2.18),  and  using  the  quasi-neutrality  approximation 

ni  “  ne  =  n  (2.19) 


we  obtain 
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-1  Ri  TbT  -i  t 


+  r  (Ri  IiJ.  +  Iel>X  £ 


(2.20) 


For  the  barium  cloud  problem  we  shall  treat  here,  we  shall  only  consider 
neutral  winds,  electric  fields,  gravity,  and  pressure  gradients  as  external 
forces.  Hence 

— ii  "  6  I±  +  mi  fij.  +  vin  mi  — nl  "  VPi/n  <2*21> 

F^i  =  -  e  Ei  +  me  gL  -  VPe/n  (2.22) 

Note  that  we  have  neglected  the  small  term  ven  rae  in  Equation  (2.22). 


We  obtain 


Ri  TbT  (e  -I.  +  mi  *j.  +  vin  mi  -Hnl  '  VPi/n> 


+RiF  fe  (1  -  Ri  )  +  Cn£  +  r^)  8l  +  Vin-iHo 


-  VP  /n  -  VP  R?  /n]  x  z 
i  e  l  — 


Since  O.Ol  <  <  1.0  we  may  neglect  ma/R^  with  respect  to  m^. 

Defining  the  Pedersen  conductivity 


_  _  „  In  nee 

aP  =  Ri  07-  W 


and  noting  that  1  -  R^1  =*  -  vln2/0i2  we  obtain 


r  mi  vp< 

J  -  a  [E  +  -i  g  +  v.  —  U  -  — - 

—I  p  L— i  e  — 1  in  e  — n.L  ne 


B  ,  _  in  _  “i“i  ffli  TT 

ml  57-  £1  +  v~i «,  +  ni  r-  5, 

1  in 


(2.23) 


(2.24) 


— (TOj  +  R-lVPJ)  x  1] 


(2.25) 


Our  need  for  an  expression  for  2j_  stems  from  our  need  for  its  divergence  to 
evaluate  V  •  J  (  -  0  by  quasi-neutrality),  as  we  shall  see  in  the  next 


section. 
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3. 


Mathematical  Model 


We  shall  model  our  physical  system  using  a  simplified  model  as 
depicted  in  Figure  2.  The  magnetic  field  lines  are  assumed  to  be  straight, 
to  be  aligned  along  the  z  axis  of  our  cartesian  coordinate  system,  and  to 
terminate  in  insulators  at  z  ■  +  «.  The  plasma  of  interest  is  threaded  by 
these  magnetic  field  lines,  and  is  divided  into  thin  planes  or  "layers”  of 
plasma  perpendicular  to  the  magnetic  field.  Since  we  have  neglected 
collisions  between  different  plasma  species,  we  may  use  the  device  of 
layers  to  treat  multiple  ion  species  at  a  single  point  in  space  simply  by 
allowing  multiple  layers  to  occupy  the  same  plane  in  space,  one  for  each 
ion  species.  In  this  way  a  "layer"  consists  only  of  a  single  ion  species 
and  its  associated  electrons. 

Our  quasi-neutrality  assumption  demands  that 


V 


J 


i? _  J  J.  J _  J 

3x  x  3y  y 


0 


(3.1) 


Integrating  Equation  (3.1)  along  z  and  noting  from  Figure  2  that  Jz 


vanishes  at  z  ■  +  «  we  obtain 


From  our  model  as  depicted  in  Figure  2  we  may  approximate  the  integral  in 
Equation  (3.2)  by  a  discrete  sum 


K 


l 

k-1 


0 


(3.4) 


where  the  subscript  k  refers  to  the  layer  number,  K  is  the  total  number  of 
layers  in  the  system,  and  Az^  is  the  thickness  of  layer  k  measured  along 
the  magnetic  field  line.  Within  a  layer,  both  the  ion  density  n  and  the 
ion-neutral  collision  frequency  u^n  are  assumed  to  be  constant  along  a 
magnetic  field  line.  This  enables  us  to  introduce  the  three  magnetic- 
field-line  integrated  quantities: 


\  5 

v 

V  s  °PkAzk  -  Vnf  RJk  W 

Ehk  H  tfl^k^pk 


(3.5) 

(3.6) 

(3.7) 


By  our  assumption  of  equipotential  magnetic  field  lines  and 
electrostatic  electric  fields 


£Ik(x»y)  ”  _  for  a11  k  (3.8) 

where  =  Eq^x  +  EQy.y  is  a  constant,  externally  imposed  electric  field. 
Then  Equation  (3.4)  becomes 


(3.9) 


where 


Text 

-Ik 


“k  •  -  fe  t\k  If)  + h  lJhk  If) 

-  -  li  3?~hk  +  ii  3?:tik 

3y  3x  +  3x  3y  (3.10) 


m.  m,  7P 
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W  W  £l  1  «“  -1 
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- -  I  7P,  +  rTVp  I 

vln  1  t  eJ 


in  „  ,  A 1 
-  - —  E  x  2, 
Cl  —a J  —  J] 


(3.11) 


and  the  subscript  k  denoting  layer  number  on  terms  within  parenthesis 
operates  on  all  terms  within  those  parentheses. 

The  system  of  equations  we  must  solve  consists  of  (3.9)— (3. 11)  and  an 
equation  of  continuity  for  each  layer: 

3Nfc 

"3t  +  V  •  [ N^v^ J  "  ®  k*  1 ,  . . . ,  K  (3.12) 


where  v^  is  given  by  Eq.  (2.6)-(2.12). 


4.  Simplified  Two-Level  Model 

We  make  the  following  further  simplifications  in  our  model: 

a)  We  find  that  for  ionospheric  barium  clouds  the  currents  parallel 

to  B  are  carried  primarily  by  electrons,  and  the  motion  of  ions  parallel  to 
B  consists  primarily  of  a  slow  diffusion  plus  a  bulk  falling  of  the  cloud 
to  lower  altitudes  (since  £  •  3  *  0).  It  is  therefore  sufficient  to 
represent  the  ions  as  an  array  of  two-dimensional  planes  or  layers  of 
plasma  perpendicular  to  B^,  each  moving  with  the  bulk  "falling”  velocity 
along  and  hence  to  treat  numerically  only  the  transport  of  ions 

perpendicular  to  within  each  layer. 

b)  We  limit  our  discussion  to  a  model  consisting  of  only  two  layers 
or  species  of  ions. 

c)  We  assume  that  the  ion  neutral  collision  frequency  is 
constant  within  a  layer. 

d)  We  assume  that  the  only  forces  acting  on  the  plasma  are  an 
external  electric  field  Eq  and  pressure  gradients. 

With  the  above  assumptions  we  have 
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nr  in— o 
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ne  v 


1  (VPi  +  Ri17PeJ]fc  x  zj 


(4.1) 


in 


If  we  assume  B  is  along  the  positive  z  axis  then  B/|B|  *  1  and  we  get 


VP, 
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(4.4) 


which  is  a  constant  within  a  layer  if  yln  is  constant  within  a  layer  (one 
of  our  assumptions).  Then 


Hpk  ■  ~  •  (constant  V(P^  +  R^P^J  x  jtj  ■  0 


(since  7^  •  (7C  x  z)  -  0  for  any  scalar  field  C). 
Our  final  equations  to  be  solved  are  then 


3N 

Tr+7i  •  (Vtu)  k->'2 


(4.5) 


(4.6) 


7i  •  (I1.!  +  to2)V]  +H-E  .  7  (r  +E  ,) 


where 


H  r  “  3x  ^Ehl  +  Zh2^3y  "  Eoy^ 
+  ly  ^Zhl  +  Eh2^H  “  Eox^ 


(4.8) 


Looking  at  Eq.  (4.7)  we  notice  that  may  be  separated  into  two  parts: 


<J>  *  4>e  +  <J>p 


(4.9) 


71  •  U'pl  +  WeI  + 


H  -  E 


•  W  +  'pi) 


(4.10) 


'.I  •  U'pl  +  Ep2)7l+pl  -  -  V'pl  -qi  +  rp2  ^ 


(4.11) 


For  reasons  which  shall  become  clear  as  we  progress,  we  shall  regard  $E  as 
that  part  of  the  potential  field  which  tends  to  drive  the  cloud  toward 
bifurcation.  In  fact,  Eq,  (4.10)  coupled  with  (4.6)  just  form  the  basic 
inviscid  two-level  equations  [Scannapieco  et  al.,  1976].  We  shall  regard 
$p  as  that  part  of  the  potential  field  which  tends  to  diffuse  or^  antl- 
dif  fuse  the  edges  of  the  cloud.  First  let  us  look  at  solutions  to 

(4.11).  If  VP^  and  VP^2  vanish  at  our  boundaries,  the  unique  solution 

consistent  with  zero  Neumann  boundary  conditions  (vanishing  of  the  normal 
derivative  of  at  the  boundary)  gives: 


^Zpl  +  Ep2^7i^P  “  Zpl  nte 


p2  n2e 


(4.12) 


In  general,  the  motion  resulting  from  the  forces  proportional  to 
7<$>D,  VP,  and  VP.,  will  be  quite  complex,  not  describable  as  the  simple 


superposition  of  a  shear  and  a  diffusion.  However,  in  the  special  case  of 
one  of  the  layers  being  uniform,  e.g.,  N0  and  Ti2  (T^  »  ion  temperature) 
hence  Z^2  an<^  ^12  uniform,  then 
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E  .  +  E  ,,  n.  e 
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(4.13) 


We  shall  now  treat  this  special  case.  One  may  conveniently  think  of  level 
1  as  representing  the  barium  cloud  and  level  2  the  uniform  background 
ionosphere. 

We  now  use  the  equation  of  state 


where  is  in  energy  units  (T^  is  the  product  of  Boltzmann's  constant  and 
the  temperature  in  degrees  Kelvin),  and  assume  thermal  conductivities  high 
enough  such  that  is  constant  and  hence 


Vpi  -  TjVn 


(4.14) 


Then 


-E  . 


T.  .Vn. 


The  velocity  in  level  1  resulting  from  the  combined  action  of  7P^  and  i p  is 
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(4.16) 


Thus 


v  ~Zp2  r.vin  c  VN, 
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Vn 

Note  we  have  used  — 
n 


VN 

N 


It  is  convenient  to  divide  v^lp  into  two  parts: 


-ilP 


ID  l 


+  -SI 


with  the  subscripts  i  and  P  suppressed  but  understood. 
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(4.19) 


That  _V£^  represents  shear  flow  at  plasma  gradients  can  be  seen  by  noting 

that  the  velocity  is  always  perpendicular  to  VN,  and  is  largest  where 
Ep2  VN 

r — -■■■-» —  ^  is  a  maximum,  decaying  in  magnitude  away  from  the  region.  The 
pi  +  ZP2  *N 

quantity  ^  is  of  interest  because  it  is  a  shear  flow  and  hence  may  be 
stabilizing  in  certain  cases  (see  Perkins  and  Doles  [1975],  Huba  et  al. 


[1982]).  We  do  not  consider  the  effects  of  _Vg  further  in  this  paper  except 

to  note  that  for  the  special  case  of  a  uniform  layer  2  being  considered 

here,  V  •  _Vg^  *  0  and  VN^  •  -  0  which  together  imply  V  •  (.N^v^J  ■  0, 

which  is  to  say  that  has  effect  on  the  time  evolution  of  Nj^.  Since 

this  breaks  the  feedback  loop,  it  is  difficult  to  see  how  a  shear 

stabilization  mechanism  could  be  active  in  this  case. 

A  consideration  of  _vp  is  really  the  heart  of  this  paper.  The  quantity 

vn  is  the  source  for  ion  diffusion  in  level  1.  More  importantly,  it  is  the 

source  for  the  diffusion  of  the  total  integrated  Pedersen  and  Hall 

conductivities  if  level  2  is  sufficiently  incompressible  (i.e., 

if  v.  /£}.  in  level  2  is  sufficiently  small),  as  we  will  show.  Let  us  look 
m  l 

at  the  effect  of  on  the  ion  continuity  equation  for  level  1: 


This  is  just  a  diffusion  equation  for  with  diffusion  coefficient 


D 

o 


(4.21) 


(It  should  be  noted  that  this  diffusion  coefficient  is  similar  to  that 
derived  by  Goldman  et  al.,  1976,  and  Vickrey  and  Kelley,  1982.) 


For  a  barium  cloud  (level  1)  at  180  km  altitude  we  take  ~  1000°  K, 
v 

— —  ~  0.06,  B  »  0.5  g,  and  get 
Ui 

Z  o 

D  -  P7  „ —  (100  m2 /sec)  (4.22) 

o  E  ,  +  Z  - 


In  the  less  dense  regions  of  the  barium  cloud  where  E  ^  is  small  compared 
to  £p2,  the  ion  diffusion  coefficient  is  quite  large  indeed.  By  contrast 
ambipolar  diffusion  rates  Induced  by  electron-ion  and  electron-neutral 
collisions  yield  diffusion  coefficients  on  the  order  of  1  m  /sec.  Thus  our 
neglect  of  these  collision  terms  seems  justified.  Unfortunately,  is  not 
the  relevant  quantity  if  one  is  interested  In  the  effect  of  _vj)  on  barium 
cloud  dynamics.  A  look  at  Eq.  (4.10)  will  convince  the  reader  that  if  one 


wants  to  affect  the  time  evolution  of  <J>E»  the  quantity  which  determines  the 


bifurcation  process,  it  is  necessary  to  change  the  quantities  T.  =  + 

Ep2  a°d  +  E^2»  c*ie  tol:al  field  line  integrated  Pedersen  and  Hall 

conductivities.  Again  we  assume  an  initially  uniform  distribution  of  Z 2 
and  N^.  We  can  only  show  the  effects  of  vq  in  the  first  instant  of  time 


here.  Once  Images  form  in  level  2,  the  problem  yields  only  to  numerical 
techniques.  We  need  the  velocity  in  level  2. 
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So  we  see  that 


a(NL  +  n2) 

at 
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(4.26) 


which  when  combined  with  quasi-neutrality,  is  simply  a  statement  that 

electrons  can't  diffuse.  This  is  consistent  with  our  neglect  of  electron- 

ion  and  electron-neutral  collisions. 

Now  let  us  look  at  the  effect  of  vn  on  E  =  E  ,  +  E  ».  A  similar 

— D  p  pi  p2 

analysis  can  be  performed  for  E  with  similar  results.  However  the 
bifurcation  effects  of  are  small  compared  to  those  of  E  ,  and  we  shall 
not  consider  them  here. 
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We  note  that  in  Eq.  (4.18)  we  may  set  VN^/N^  =  ^pl/^pl  so  *-^at 
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Now  since 
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we  have 
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2 

nominal  100  m  vsec  diffusion  coefficient  is  not  only  multiplied  by 

+  £p2J»  re<*ucing  its  effectiveness  in  the  dense  portions  of  the 
cloud,  but  it  is  further  multiplied  by  the  factor  F  where 


vinRi/nlh 

(VWi 


(4.37) 


Note  that  if  v.  R, /n.  is  the  same  for  both  layers,  the  diffusion  is  reduced 
in  l  l 

to  zero.  Furthermore,  if  1^2  ^  ^VinRi^i^l  compressible 
background  plasma,  e.g.,  an  E  region)  the  diffusion  coefficient  is 
negative .  That  is,  gradients  in  will  actually  be  steepened,  rather  than 
smoothed.  For  this  reason  we  conclude  that  unless  some  rather  fast 
recombination  chemistry  is  taking  place  in  the  E  region,  the  presence  of  a 
background  E  region  is  destabilizing,  in  accord  with  the  findings  of 
Francis  and  Perkins  [1975]. 

The  opposite  extreme  of  J  2  (.  v^nRj_/ft  jj  (background  F  region 
ionosphere)  is  more  interesting  since  it  closely  approximates  the 
conditions  of  most  barium  releases.  It  also  yields  the  full  nominal  100 
m  /sec  diffusion  coefficient,  since  F  a  1.0,  giving  us  some  hope  of 


stopping  the  bifurcation  process  at  scale  sizes  of  hundreds  of  meters. 

Note  that  in  this  case  the  diffusion  rate  for  Z  is  the  same  as  that  for 

P 

N^.  Also  note  that  the  case  of  10 1^2  sma^  is  t*ie  case  where  the 
Pedersen  mobility  and  hence  compressibility  of  the  background  is  extremely 
low,  meaning  that  only  minimal  images  will  be  formed.  It  is  a  good 

approximation  to  take  to  be  uniform  for  all  time.  One  might  well  ask 
how  this  is  possible  given  that  the  electrons  can't  diffuse,  and  that  the 
barium  ions  are  being  allowed  to  diffuse  in  level  1.  The  answer  is  that 
there  are  (small)  images:  the  electrons  lost  to  a  field  line  in  level  1 

are  "soaked  up"  on  the  same  field  line  in  level  2.  But  the  number  of 
electrons  N2  in  level  2  is  extremely  large.  (N2  must  be  large  in  order  to 
contribute  a  significant  Pedersen  conductivity  and  yet  have 
small  (see  Eq.  2.24)).  Thus  the  electrons  that  are  soaked  up  in 

level  2  induce  only  an  extremely  small  percentage  change  in  Z 


5.  Diffusion  Characteristics  for  a  Simplified  Two-level  Model  with  a 
Nearly  Incompressible  Background  Layer 

In  the  last  section  we  demonstrated  that  when  the  background 

ionosphere  (layer  2)  is  nearly  incompressible  C ^ 1^2  is  sma1^)  then 

the  primary  effect  of  a  finite  barium  ion  temperature  is  to  introduce  into 

the  continuity  equation  for  both  N.  and  for  Z  a  diffusion  term  with 

1  P 

diffusion  coefficient  DQ  given  by  Eq.  (4.21).  For  parameters  appropriate 
to  a  barium  cloud  at  180  km  altitude  we  get  Eq.  (4.22): 
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Since  level  2  i3  nearly  incompressible, 
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we  may  regard  T. 


as  a  constant. 


a 


Since  it  is  proportional  to  N^,  the  integrated  barium  cloud  density,  E  ^ 

is  larger  in  the  center  of  the  cloud  and  decays  toward  zero  at  the  cloud 

edges.  This  means  that  the  diffusion  coefficient  approaches  100  m~/sec 

near  the  cloud  edges,  but  may  be  considerably  less  than  that  value  near  the 

cloud  center  if  Z  .  »  E  -• 
pi  p2 

Previous  attempts  to  explain  the  persistence  of  scale  sizes  of 
hundreds  of  meters  for  long  periods  of  time  after  barium  releases  have 
inevitably  encountered  the  following  problems:  1)  The  diffusion 

coefficients  required  to  stop  the  bifurcation  of  plasma  clouds  several 
hundred  meters  in  diameter  were  higher  than  one  could  explain 
theoretically;  2)  assuming  the  existence  of  such  large  (ordinary)  diffusion 
coefficients,  their  effect  would  not  only  be  to  stop  the  bifurcation  of  the 
cloud,  but  also  to  diffuse  the  cloud  away  entirely  in  several  minutes.  The 
form  of  the  diffusion  coefficient  given  by  Eq.  (4.22)  would  seem  to  offer 
us  some  hope  of  overcoming  both  of  these  problems.  Its  success  depends  on 
the  accuracy  of  the  following  qualitative  picture  of  barium  cloud 
bifurcation:  We  postulate  a  barium  cloud  such  that  E^  >  5  E^  so  that  the 
diffusion  well  inside  the  cloud  is  considerably  reduced  over  that  at  its 
edges.  We  further  postulate  that  bifurcation  is  a  process  which,  in  the 
absence  of  diffusion,  starts  at  the  edges  of  barium  clouds  and  works  its 
way  inward.  Our  diffusion  then  satisfies  the  needed  criteria  perfectly: 
1)  It  is  very  effective  at  the  edges  of  the  cloud,  stopping  bifurcation 
before  it  starts;  2)  it  is  very  ineffective  at  the  central  core  of  the 
cloud,  giving  it  a  long  lifetime.  Note  that  if  this  picture  is  correct 
then  the  barium  clouds  commonly  known  as  large  M  clouds  where 
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will  have  extremely  long  lifetimes.  In  fact  for  M  >  100,  the  decay  of  the 
core  is  probably  determined  by  ambipolar  diffusion  rates.  The  accuracy  of 
this  picture  is  best  determined  by  numerical  simulation  techniques,  to 
which  we  turn  in  the  next  section. 

6.  Numerical  Simulations  Using  a  Simplified  Two-level  Model  with  a  Nearly 

Incompressible  Background  Layer 

In  this  section  we  attempt  to  answer  numerically  the  question  posed  in 
the  last  section:  if  we  consider  a  barium  cloud  whose  background 
conductivity  is  due  primarily  to  high  F  region  plasma  (nearly 
incompressible  plasma),  is  the  diffusion  coefficient  given  by  Eq.  (4.22) 
sufficient  to  stop  the  bifurcation  of  that  cloud  at  a  diameter  of  several 
hundred  meters?  In  order  to  isolate  the  effects  of  this  diffusion 
coefficient,  we  will  make  some  further  simplifications  in  our  two-level 
model:  1)  We  assume  that  N2  and  hence  remain  uniform  during  all  times 
of  interest.  This  is  completely  consistent  with  our  assumption  of  a  nearly 
incompressible  background  plasma,  as  discussed  at  the  end  of  Section  4;  2) 
We  neglect  the  shear  component  v^  of  the  pressure-induced  ion  velocity  as 
defined  in  Eq.  (4.19).  Again  this  is  consistent  with  our  assumption  of  a 
nearly  incompressible  background  plasma  which  in  turn  implies  a 
uniform  E^*  33  shown  in  the  discussion  following  Eq.  (4.19);  3)  VJe 
neglect  the  "Hall  terra"  H  in  Eq .  (4.4)  and  hence  in  Eq.  (4.10).  This  is  a 
good  approximation  as  long  as  v ^n/f2 ^  <0.1;  4)  We  neglect  Pedersen 
convection.  That  is,  for  all  terms  except  those  causing  ion  diffusion,  we 
approximate  the  ion  velocity  with  the  electron  velocity  (see  Eq.  (6.4) 


below).  We  know  fron  our  previous  study  [Zalesak  et .  al , ,  1983]  that  the 


Pedersen  mobility  of  the  ions  in  response  to  (E^  -  V<fc E )  can  have  an  effect 
on  the  "freezing"  phenomenom,  and  it  is  our  desire  to  isolate  the  effects 
of  the  pressure-induced  diffusion.  Hence  we  neglect  ion  Pedersen  mobility 


except  in  the  diffusion  term. 

The  final  equations  to  be  solved  numerically  are  then 

3N 

yp  +  V)-(N1v1)  -  V[.(D7N1) 

(6.1) 

3N_ 

2  -  0 

3t 

(6.2) 

V<Epl  +  V  W  -  I  o*W  +  V 

(6.3) 

h  m  W  (Io  '  7V  x  1 

(6.4) 

D  -  D 

E  ,  +  E  _  io 
pi  p2 

(6.5) 

„  _  »T  f  Vin^i  ,  ce 

EPk  =  VTTTvln7a1)^k  W 

(6.6) 

The  initial  conditions  for  the  simulations  were  as  follows 

Zp2  “  1*0  (6.7) 

E  .  -  4.0  exp  (-r4/ru)  (6.8) 

pi  o 

r2  **  (x  -  x  )2  +  (y  -  y  )2  (6.9) 


rQ  *  (250  m)  (1  +  0.001  sin  80) 


(6.10) 


9  •  arctan  ((y-yQ)/(x  -  xq)J 


(6.11) 


B  ■  0.5  gauss 


(6.12) 


cE  /B  -  100  m/sec 
o 


(6.13) 


where  xQ  and  yQ  are  the  coordinates  of  the  cloud  center  and 
the  x  B  velocity  (cEQ/B)  is  in  the  negative  y  direction. 

The  above  equations  are  solved  numerically  on  a  finite  difference  grid 
in  an  x-y  cartesian  geometry  perpendicular  to  the  magnetic  field,  which  is 
assumed  to  be  aligned  along  the  z-axis.  The  83  x  60  grid  is  stretched  in 
all  four  directions,  with  the  central  56  x  32  portion  of  the  grid,  which  is 
centered  on  the  steepening  backside  of  the  cloud,  having  a  grid  spacing  of 
10  m  in  both  directions.  The  grid  stretching  allows  the  boundaries  to  be 
placed  A  km  away  from  the  edges  of  the  central  uniform  mesh  in  all 
directions.  The  potential  equation  (6.3)  is  solved  using  a  vectorized 
incomplete  Cholesky  conjugate  gradient  algorithm  due  to  Hain  (1980),  which 
is  an  extension  of  the  algorithm  of  Kershaw  (1978).  The  continuity 
equation  (6.1)  is  integrated  forward  in  time  using  the  fully 
multidimensional  flux-corrected  transport  algorithms  of  Zalesak  [1979]. 


In  Figure  3  we  show  isodensity  contours  of  E  for  the  initial 


conditions  for  all  of  the  calculations  we  shall  show.  We  have  performed 


three  calculations,  varying  the  value  of  from  zero  up  to  the  physically 

2 

realistic  value  of  100  m  /sec.  The  purpose  of  this  sequence  is  to  show 


first  that  the  diffusion  does  have  an  effect  on  the  evolution  of  the  cloud, 


and  second  to  show  that  this  effect  is  much  larger  than  that  due  to  any 

numerical  diffusion  which  may  be  present  in  the  calculations.  In  Figure  4 

we  show  the  results  for  D,  *  0  .  Note  that  even  at  the  very  early  time  of 

14  seconds,  bifurcation  has  already  started.  By  24  seconds,  the  main 

portion  of  the  cloud  has  completely  sheared  in  two.  In  Figure  5,  we  show 

2 

the  results  for  D  *  25  m  /sec  ,  one  fourth  of  the  physically  realistic 

value.  The  development  of  the  cloud  has  slowed  considerably,  showing  that 

even  at  one  fourth  its  physically  realistic  value,  the  physical  diffusion 

is  dominating  the  effects  of  whatever  numerical  diffusion  may  be  present. 

At  14  seconds  the  cloud  has  not  clearly  begun  to  bifurcate.  At  25  seconds, 

bifurcation  has  occurred,  but  is  qualitatively  and  quantitatively  different 

from  that  which  took  place  when  D^Q  =*0.  In  Figure  6  we  show  the  evolution 

of  the  cloud  with  the  full,  physically  realistic  value  for  Dio  of  100 

m  /sec.  Bifurcation  has  now  been  halted.  Rather,  the  cloud  evolves  into  a 

streamlined  "bullet"  or  "tadpole"  shape  in  time,  and  then  undergoes  a  slow 

diffusive  decay.  Even  at  the  latest  time  calculated,  75  sec,  there  is  no 

hint  whatsoever  of  imminent  bifurcation.  Noting  that  for  this  cloud 

M  ■  (E  ,  +  E  „)/E  -  is  5  In  the  center  of  the  cloud,  which  yields  an 
pi  pZ  pZ 

effective  diffusion  coefficient  of  20  m4/sec  there.  This  is  still  enough 
to  result  In  substantial  decay  of  the  central  core  on  time  scales  of 
minutes.  Clouds  with  larger  values  of  M  in  their  centers  would  of  course 
be  longer-lived. 


7. 


Conclusions  and  Future  Work 


The  primary  conclusion  of  this  paper  can  be  stated  simply:  Under 

realistic  ionospheric  conditions  (where  the  plasma  cloud  couples  to  an 
incompressible  background  i.e.,  in  F  region),  barium  striations  hundreds  of 
meters  in  diameter  can  be  long-lived,  quasi-stable,  non-bifurcating 
structures.  These  structures  may  resemble  "tadpoles",  with  a  dense  head, 
steep  density  gradients  at  the  front,  and  a  long,  less  dense  tail.  If  what 
one  means  by  the  term  "freezing"  is  the  above  phenomenon,  then  we  have 
shown  that  freezing  does  indeed  exist. 

The  next  obvious  step  is  to  test  some  of  the  rpproximations  we  have 
made  in  our  preliminary  numerical  simulation  model.  Pedersen  convection 
should  be  put  back  into  the  model,  as  well  as  the  "Hall  terms”  in  the 
potential  equation.  Finite  images  should  be  accounted  for.  Also  we  would 
like  to  test  the  combined  effect  of  the  presence  of  both  a  relatively 
incompressible  F  region  background  and  a  compressible  E  region  on  the 
dynamics  of  the  cloud.  Recombination  chemistry  must  be  included  in  this  E 
region  if  it  is  to  be  modeled  accurately  (see  Vickery  and  Kelley,  1982]. 

Another  of  our  approximations  which  should  be  studied  in  detail  is 
that  the  barium  cloud  and  background  ionosphere  are  being  driven  by  an 
externally  imposed  electric  field,  or  equivalently,  by  a  neutral  wind  which 
doesn't  vary  in  altitude.  Were  we  to  address  an  altitude-dependent  neutral 
wind  as  a  driver,  some  of  our  simple  results  would  become  less  so.  The 
most  important  factor  would  be  that  one  could  no  longer  argue  that  the 
bifurcation  of  the  cloud  depended  only  on  and  since  the  potential 

equation  (4.10)  would  then  contain  terms  which  depended  explicitly  on 
and  N?.  Although  this  would  not  change  the  results  for  the  incompressible 
background  ionosphere  considered  here,  it  would  appear  to  allow  for  the 
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possibility  of  this  end-shorting-induced  diffusion  being  an  effective 
bifurcation  inhibiting  mechanistr  even  when  the  background  ionosphere  is 
compressible.  We  shall  consider  this  case  in  our  future  work  in  tnis  area. 

It  is  also  clear  that  we  must  address  the  validity  of  the  perfect 
mapping  of  the  electrostatic  potential  along  B.  Given  that  we  are  now 
dealing  with  structures  only  a  few  hundred  meters  in  diameter,  with 
gradient  scale  lengths  of  a  few  tens  of  meters,  it  seems  likely  that  we 
will  have  to  deal  with  the  question  of  how  far  the  polarization  field 
produced  by  the  barium  cloud  maps  along  the  magnetic  field  (Goldman  et.  al, 
1976,  Glassman  and  Sperling,  1983).  One  possible  effect  might  be  that  the 
effective  value  of  M  would  be  increased,  since  the  barium  cloud  could  no 
longer  "see"  background  plasma  that  was  far  from  the  cloud  position  along 
B. 
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Figure  1.  Sketch  of  the  time  evolution  of  a  typical  barium  cloud  in  a 
plane  perpendicular  to  the  magnetic  field,  subject  to  an  upward 
directed  neutral  wind  or  equivalently  to  a  rightward  directed 
external  electric  field.  Lines  demarking  the  cloud  denote 
plasma  density  contours,  with  the  highest  plasma  density  in  the 


center  of  cloud. 


PLASMA 


Figure  2.  Model  of  plasma  and  magnetic  field  geometry  used  in  this 
paper.  Field  lines  terminate  on  insulators  at  z  =■  +®.  Plasma 
is  divided  into  "layers”  along  z  for  mathematical  and  numerical 
treatment.  Each  layer  consists  of  a  single  ion  species  and  its 
associated  electrons.  Multiple  collocated  ion  species  are 
treated  by  having  multiple  collocated  "layers". 
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Figure  3 


Isodensity  contours,  with  areas  between  alternate  contours  shown 
cross-hatched,  for  the  barium  cloud  used  for  the  three  numerical 
simulations  shown  in  this  paper.  The  center  of  the  cloud  is 
located  at  the  point  (0,0).  The  boundaries  in  x  are  located 
at  +  4.9  km.  The  boundaries  in  y  are  located  at  -4.0  km  and 
+5.3  km  respectively.  The  cross-hatching  is  done  at  every  other 
grid  line,  so  that  one  rectangle  of  cross  hatching  represents  a 
2x2  array  of  computational  cells.  Contours  for  this  and  all 


other  plots  in  this  paper  are  drawn  at  values  for  Z  ^  of  0.31, 


0.71,  1.24,  1.92,  and  2.82  (logarithmic  spacing  for  total 

Pedersen  conductivity). 
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Figure  5.  Time  evolution  of  the  cloud  depicted  In  Figure  3  for  the 

2 

case  =  25  m  /sec  .  Shown  are  lsodensity  contours  at  a)  14 
sec,  b)  24  sec.  Contour  values  are  as  given  in  Figure  3. 
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Figure  6.  Time  evolution  of  the  cloud  depicted  in  Figure  3 

2 

case  D.  =  100  m  /sec  (the  physically  realistic  case), 
io 

are  isodensity  contours  at  a)  14  sec,  b)  25  sec,  c)  39 
50  sec,  e)  64  sec,  f)  75  sec. 
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Figure  6,  cont'd.  Time  evolution  of  the  cloud  depicted  in  Figure  3  for 

2 

the  case  *  100  m  /sec  (the  physically  realistic  cast 
Shown  are  isodensity  contours  at  a)  14  sec,  b)  25  sec,  c) 
sec,  d)  50  sec,  e)  64  sec,  f)  75  sec. 


Figure  6,  cont'd. 


Time  evolution  of  the  cloud  depicted  in  Figure  3  for 
2 

the  case  D,  ■  100  m  /sec  (the  physically  realistic  case 
lo 

Shown  are  isodensity  contours  at  a)  14  sec,  b)  25  sec,  c) 
sec,  d)  50  sec,  e)  64  sec,  f)  75  sec. 
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COMMANDER 

U.S.  ARMY  COMM-ELEC  ENGRG  INSTAL  AGY 
FT.  HUACHUCA,  AZ  85613 

01CY  ATTN  CCC-EMEO  GEORGE  LANE 

COMMANDER 

U.S.  ARMY  FOREIGN  SCIENCE  &  TECH  CTR 
220  7TH  STREET,  NE 
CHARLOTTESVILLE,  VA  22901 
01CY  ATTN  DRXST-SD 

COMMANDER 

U.S.  ARMY  MATERIAL  DEV  &  READINESS  CMD 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

01CY  ATTN  DRCLDC  J.A.  BENDER 

COMMANDER 

U.S.  ARMY  NUCLEAR  AND  CHEMICAL  AGENCY 
7500  BACKLICK  ROAD 
BLDG  2073 

SPRINGFIELD,  VA  22150 
01CY  ATTN  LIBRARY 


U.S.  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND,  MD  21005 

01CY  ATTN  TECH  LIBRARY  EDWARD  3AICY 

COMMANDER 

U.S.  ARMY  SATCOM  AGENCY 
FT.  MONMOUTH,  NJ  07703 

01CY  ATTN  DOCUMENT  CONTROL 

COMMANDER 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  ARSENAL,  AL  35809 
01CY  ATTN  JIM  GAMBLE 

DIRECTOR 

U.S.  ARMY  TRADOC  SYSTEMS  ANALYSIS  ACTIVITY 
WHITE  SANDS  MISSILE  RANGE,  NM  88002 
01CY  ATTN  ATAA-SA 
01CY  ATTN  TCC/F .  PAYAN  JR. 

01CY  ATTN  ATTA-TAC  LTC  J.  HESSE 

COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
WASHINGTON,  D.C.  20360 

01CY  ATTN  NAVALEX  034  T.  HUGHES 
01CY  ATTN  PME  117 
01CY  ATTN  PME  117-T 
01CY  ATTN  CODE  5Cil 

COMMANDING  OFFICER 
NAVAL  INTELLIGENCE  SUPPORT  CTR 
4301  SUITLAND  ROAD,  BLDG.  5 
WASHINGTON,  D.C.  20390 

01CY  ATTN  MR.  DUBBIN  STIC  12 

01CY  ATTN  NISC-50 

01CY  ATTN  CODE  5404  J.  GALET 

COMMANDER 

NAVAL  OCCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 
01CY  ATTN  J.  FERGUSON 
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NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 


01CY 

ATTN 

CODE 

4700 

S. 

L.  Ossakow 

26  CYS  IF 

UNCLASS. 

1  CY  IF  CLASS) 

01CY 

ATTN 

CODE 

4701 

I 

Vitkovltsky 

01CY 

ATTN 

CODE 

4780 

J. 

Huba  (100 

cys  : 

IF  UNCLASS, 

1 

CY  IF  CLASS) 

01CY 

ATTN 

CODE 

7500 

01CY 

ATTN 

CODE 

7550 

01CY 

ATTN 

CODE 

7580 

01CY 

ATTN 

CODE 

7551 

01CY 

ATTN 

CODE 

7555 

01CY 

ATTN 

CODE 

4730 

E. 

MCLEAN 

01CY 

ATTN 

CODE 

4108 

01CY 

ATTN 

CODE 

4730 

B. 

RIPIN 

20CY 

ATTN 

CODE 

2628 

COMMANDER 

NAVAL  SEA  SYSTEMS  COMMAND 
WASHINGTON,  D.C.  20362 
01CY  ATTN  CAPT  R.  PITKIN 

COMMANDER 

NAVAL  SPACE  SURVEILLANCE  SYSTEM 
DAHLGREN,  VA  22448 

01CY  ATTN  CAPT  J.H.  BURTON 

OFFICER-IN-CHARGE 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING,  MD  20910 
01CY  ATTN  CODE  F31 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20376 
01CY  ATTN  NSP-2141 
01CY  ATTN  NSSP-2722  FRED  WIMBERLY 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
DAHLGREN,  VA  22448 

01CY  ATTN  CODE  DF-14  R.  BUTLER 

OFFICER  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 
01 CY  ATTN  CODE  465 


COMMANDER 

AEROSPACE  DEFENSE  COMMAND/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB ,  CO  80912 

01CY  ATTN  DC  MR.  LONG 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND/ XPD 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 

01CY  ATTN  XPDQQ 

01CY  ATTN  XP 

AIR  FORCE  GEOPHYSICS  LABORATORY 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  OPR  HAROLD  GARDNER 
01CY  ATTN  LKB  KENNETH  S.W.  CHAMPION 
01CY  ATTN  OPR  ALVA  T.  STAIR 

01CY  ATTN  PHD  JURGEN  BUCHAU 

01CY  ATTN  PHD  JOHN  P.  MULLEN 

AF  WEAPONS  LABORATORY 
KIRTLAND  AFT,  NM  87117 
01CY  ATTN  SUL 

01CY  ATTN  CA  ARTHUR  H.  GUENTHER 
01CY  ATTN  NTYCE  1LT.  G.  KRAJEI 

AFTaC 

PATRICK  AFB,  FL  32925 
01CY  ATTN  TF/MAJ  WILEY 
01CY  ATTN  TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT -PATTERSON  AFB,  OH  45433 
01CY  ATTN  AAD  WADE  HUNT 
01CY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  6  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20330 
01CY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 
01CY  ATTN  J.  DEAS 


01CY 

ATTN 

CODE 

461 

01CY 

ATTN 

CODE 

402 

HEADQUARTERS 

% 

OICY 

ATTN 

CODE 

420 

ELECTRONIC  SYSTEMS  DI  VISION/ YSEA 

01CY 

ATTN 

CODE 

421 

DEPARTMENT  OF  THE  ATR  FORCE 

HANSCOM  AFB,  MA  01732 
01CY  ATTN  YSEA 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AF8,  MA  01731 

01CY  ATTN  DCKC  MAJ  J.C.  CLARK 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  NICD  LIBRARY 

01CY  ATTN  ETDP  B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
CRIFFISS  AFB,  NY  13441 

01CY  ATTN  DOC  LIBRARY/TSLD 
01CY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZJ 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

01CY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
01CY  ATTN  NRT 

01CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01 CY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

01CY  ATTN  MNNL 

COMMANDER 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

01CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

01CY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 

01CY  ATTN  DOC  CON  FOR  L-389  R.  OTT 

01CY  ATTN  DOC  CON  FOR  L-31  R.  RACER 

01CY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 

01CY  ATTN  DOC  CON  FOR  J.  WOLCOTT 
01CY  ATTN  DOC  CON  FOR  R.F.  TASCHEK 


01CY 

ATTN 

DOC 

CON 

FOR 

E. 

JONES 

01CY 

ATTN 

DOC 

CON 

FOR 

J. 

MALIK 

01CY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

01CY 

ATTN 

DOC 

CON 

FOR 

J. 

ZINN 

01CY 

ATTN 

DOC 

CON 

FOR 

P. 

KEATON 

01CY 

ATTN 

DOC 

CON 

FOR 

D. 

WESTER  VELT 

01CY 

ATTN 

D.  SAPPENFIELD 

SANDIA 

LABORATORIES 

P.O.  BOX  5800 

ALBUQUERQUE, 

NM 

87115 

01CY 

ATTN 

DOC 

CON 

FOR 

W. 

BROWN 

01CY 

ATTN 

DOC 

CON 

FOR 

A. 

THORNBROUCH 

01CY 

ATTN 

DOC 

CON 

FOR 

T. 

WRIGHT 

01CY 

ATTN 

DOC 

CON 

FOR 

D. 

DAHLGREN 

01CY 

ATTN 

DOC 

CON 

FOR 

3141 

01CY 

ATTN 

DOC 

CON 

FOR 

SPACE  PROJECT  1 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC  SANDIA  LABORATORIES 

HANSCOM  AFB,  MA  01731  LIVERMORE  LABORATORY 

01CY  ATTN  EEP  A.  LORENTZEN  P.O.  BOX  969 

LIVERMORE,  CA  94550 

DEPARTMENT  OF  ENERGY  01CY  ATTN  DOC  CON  FOR  B.  MURPHEY 

LIBRARY  ROOM  G-042  01CY  ATTN  DOC  CON  FOR  T.  COOK 

WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  FOR  A.  LABOWITZ  OFFICE  OF  MILITARY  APPLICATION 

DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  DR.  YO  SONG 
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OTHER  GOVERNMENT 


DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

OICY  (ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

OICY  ATTN  W.  UTLAUT 
01CY  ATTN  D.  CROMBIE 
OICY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  4  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
OICY  ATTN  R.  GRUBB 
OICY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
OICY  ATTN  I.  GARFUNKEL 
OICY  ATTN  T.  SALMI 
OICY  ATTN  V.  JOSEPHSON 
OICY  ATTN  S.  BOWER 
OICY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

OICY  ATTN  L.  SLOAN 
OICY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
OICY  ATTN  C.  PRETTIE 
OICY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 


SEATTLE 

,  WA 

98124 

OICY 

ATTN 

G.  KEISTER 

OICY 

ATTN 

D.  MURRAY 

OICY 

ATTN 

G.  HALL 

OICY 

ATTN 

J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  D.B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOGSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wilshire  Blvd. 

Santa  Monica,  Calif  90401 
OICY  ATTN  C.B.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
OICY  ATTN  F.  REIBERT 
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GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 
H.MES 

COURT  STREET 
SYRACUSE,  NY  13201 

OICY  ATTN  G.  MILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
OICY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
OICY  ATTN  TECHNICAL  LIBRARY 
OICY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  CRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  DICK  STEINHOF 

HSS ,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
130  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
OICY  ATTN  J.M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFARD 
OICY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

OICY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DAS  I  AC 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  4  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  6  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  J.B.  CLAD IS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 
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MCD3NNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R.W.  HALPRIN 
01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
7  35  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CY  ATTN  P.  FISCHER 
01CY  ATTN  W.F.  CREVTER 
01CY  ATTN  STEVEN  L.  GUTSCHE 
01CY  ATTN  R.  BOGUSCH 
01CY  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KILB 
01CY  ATTN  DAVE  SOWLE 
01CY  ATTN  F.  FA JEN 
01CY  ATTN  M.  SCHEIBE 
OICY  ATTN  CONRAD  L.  LONGMIRE 
01CY  ATTN  B.  WHITE 

MISSION  RESEARCH  COP.P. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NEW  MEXICO  87106 
OICY  R.  STELLINCWERF 
OICY  M.  ALME 
OICY  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

OICY  ATTN  JOHN  MORGANSTERN 
OICY  ATTN  G.  HARDING 
OICY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

OICY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 


OICY 

ATTN 

FORREST  GILMORE 

OICY 

ATTN 

WILLIAM  B.  WRIGHT, 

OICY 

ATTN 

ROBERT  F.  LELEVIER 

OICY 

ATTN 

WILLIAM  J.  KARZAS 

OICY 

ATTN 

H.  ORY 

OICY 

ATTN 

C.  MACDONALD 

OICY 

ATTN 

R.  TURCO 

OICY 

ATTN 

L.  DeRAND 

OICY 

ATTN 

W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

OICY  ATTN  VINCE  TRAPANI 
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SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

01CY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVEN  SNOOD  AVENUE 
MENLO  PARK,  CA  94025 


OICY 

ATTN 

DONALD  NEILSON 

OICY 

ATTN 

ALAN  BURNS 

OICY 

ATTN 

G.  SMITH 

OICY 

ATTN 

R.  TSUNODA 

OICY 

ATTN 

DAVID  A.  JOHNSON 

OICY 

ATTN 

WALTER  G.  CHESNUT 

OICY 

ATTN 

CHARLES  L.  RINO 

OICY 

ATTN 

WALTER  JAYE 

OICY 

ATTN 

J.  VICKREY 

OICY 

ATTN 

RAY  L.  LEADA3RAND 

OICY 

ATTN 

C.  CARPENTER 

OICY 

ATTN 

G.  PRICE 

OICY 

ATTN 

R.  LIVINGSTON 

OICY 

ATTN 

V.  GONZALES 

OICY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,  JR. 
01CY  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  6  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 
OICY  ATTN  D /  STOCKWELL 
SNTF/1575 

VT  SI  DYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  NASS  01803 
OICY  ATTN  VI.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 


END 

FILMED 

9-84 


DTIC 


